Fingerprints of Spin-Orbital Physics in Crystalline O2 
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The alkali hyperoxide KO2 is a molecular analog of strongly-correlated systems, comprising of 
orbitally degenerate magnetic OJT ions. Using first-principles electronic structure calculations, we 
set up an effective spin-orbital model for the low-energy molecular orbitals and argue that many 
anomalous properties of KO2 replicate the status of its orbital system in various temperature regimes. 

PACS numbers: 75.50.Xx, 71.10.Fd, 75.30.Et, 71.15.-m 



Magnetic substances without d- or /-elements are 
exotic. Several rare examples include sp-impurities 
in alkali-metal hosts [H, CaBg Q, nonstoicheometric 
CaO as well as some carbon-based materials Be- 
sides them, the solid oxygen is an absolutely unique mag- 
netic system. This is because O2 is the only elementary 
molecule, whose ground state is different form the con- 
ventional spin singlet. Therefore, if these molecules form 
a crystal (either by cooling or by pressurizing), it may be- 
come magnetic. In fact, the antiferromagnetism of solid 
O2 is well known @, @, 0| , and can be anticipated from the 
electronic structure of the single O2 molecule, having two 
unpaired electrons in a doubly-degenerate n g shell |6j. 

Since O2 is a good oxidizer and can easily take an addi- 
tional electron when it brought in contact with alkali ele- 
ments, there is another way of making the crystalline O2, 
in a form of ionic crystals (called "alkali hyperoxides" ) , 
which arc similar to NaCl, but distorted due to elastic 
interactions involving O2 molecules. One typical exam- 
ple is KO2, which around the room temperature forms 
the body-centered tetragonal (bet) lattice (Fig. [1]) [1, 
The uniqueness of this situation is that O2 has an extra 
electron residing in the doubly-degenerate n g shell, and 
KO2 is the first molecular crystal, where the magnetism 
not only exists, but interplays with the orbital degrees of 
freedom, like in famous KCuF 3 andLaMn0 3 [3,[iH- 
Typically, such materials are known for their rich mag- 
netic properties, which depend on the orbital state. 
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FIG. 1: (Color online) Tetragonal phase of KO2 (left) and 
electronic structure in LDA (right). Left inset shows results 
of spin-polarized calculations. Right inset shows the 7r g bands 
near the Fermi level (located at zero energy). 



This is clearly the case also for K02- The Curie- Weiss 
temperature (Tew), extracted from the behavior of mag- 
netic susceptibility (xm) in two paramagnetic regions 
changes the sign at around 200 K [f|, indicating at the 
change of character of intermolecular interactions, from 
antifcrromagnetic (AFM) to ferromagnetic (FM) with de- 
crease of the temperature (T). The role of the crystal 
distortion in this region is somewhat controversial. On 
the one hand, there are some indications that the actual 
symmetry is lower than bet [1, 0] . On the other hand, 
no profound change of the magnetic behavior has been 
observed. For example, a sizable orbital contribution to 
the magnetization persists down to 13 K. Therefore, even 
if exists, the crystal distortion does not seem to fully 
control the magnetic behavior of KO2 in the paramag- 
netic region. If so, what causes the change of Tew 1 - An- 
other interesting point is that below 13 K the situation is 
completely different and the behavior of KO2 is strongly 
affected by frozen-in rotations of O2 molecules, which 
quenche the orbital magnetization and stabilize an AFM 
order. Nevertheless, the transition temperature is sur- 
prisingly low. Moreover, the AFM order can be broken 
up by a magnetic field. The transition is accompanied 
by structural changes, and the effect is called "magneto- 
gyration" [121 ]. 

The purpose of this work is to establish a theoretical 
basis for understanding the fascinating magnetic prop- 
erties of KO2 and other alkali hyperoxides. Using first- 
principles electronic structure calculations, we set up an 
effective spin-orbital model for the low-energy molecular 
states. Then, we show that the magnetic phase diagram 
of KO2 replicates the status of its orbital system, which 
evolves (under heating) from a quenched state, to a re- 
gion of relativistic spin-orbital correlations, and finally to 
the region of independent spin and orbital disorder. 

The O2 molecule appears to be the building block not 
only of the crystal, but also of the electronic structure 
of KO2. The hybridization within a single molecule is 
so strong that it leads to the splitting and formation of 
quite distinguishable molecular levels. The interaction 
between molecules is considerably weaker, so that the 
molecular orbitals form a group of narrow nonoverlapping 
bands (Fig. [TJ we use the experimental lattice parame- 
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ters reported in Ref. [§]). Thus, there is a clear analogy 
with atomic limit in the theory of strongly-correlated sys- 
tems [T3| , except that now the molecular orbitals play the 
same role as atomic orbitals in a conventional case. 

According to electronic structure calculations in the 
local- (spin)-density approximation [L(S)DA], the doubly- 
degenerate "Kg band located near the Fermi level is formed 
by antibonding molecular p x and p y orbitals (Fig. [1]) [l3j . 
The TT g bandwidth is comparable with the exchange spin- 
splitting (A ex ~l eV). Hence, the system is half-metallic, 
and its magnetic properties can be understood in terms of 
FM double exchange (D) and AFM superexchange (S) 
interactions, which can be expressed through the first 
and second moments of occupied density of states for the 



1-spin TT g band as J 



/2z and J s =-m^/zA c 



where z is the coordination number 11[ . Then, the Curie 



temperature will also consist of two parts, Tc=T^+T^,, 
which can be evaluated in the mean-field approxima- 
tion Q. This yields Tg=-Am^ /15k B « 882 K and 
T ( f=-m( 2 )/3fc s A ox « -558 K. Thus, T c « 324 K, and 
KO2 is expected to be a good ferromagnet, being in 
straight contrast with the experimental data. The root 
of the problem is the Coulomb correlations in the narrow 
TTg band, which are greatly oversimplified in LSDA. 

In order to treat these effects rigorously, we derive an 
effective Hubbard-type model for the ir g band, starting 
from the LDA band structure: 
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The method has been explained in details in Ref. [1 51 ] . 
The basic difference here is that each lattice-point (i or 
j) corresponds to an O2 molecule (rather than to a single 
atomic site). Each Greek symbol stands for a combina- 
tion (to, s) of orbital (m= p x or p y ) and spin (s= f or 
I) indices. Each orbital index refers to the antibonding 
molecular p- orbital of either x (m=l) or y (to=2) symme- 
try. The off-diagonal (with respect to i and j) elements 
of hff stand for transfer integrals, which do not depend 



on spin indices: i.e., ft^ = ||/i£r || = ||i™ m \\6 SS > for any com- 
bination of a=(m,s) and [3=(m' , a ). They are derived 
from the LDA band structure using an exact version of 
the downfolding method [HI, [l6| • The results are sum- 
marized in Table |T1 up to fourth nearest neighbor. Other 
parameters are considerably smaller. The diagonal ma- 
trix elements, =(cv|£L • S|/3), stand for the spin-orbit 
interaction (SOI) is a single molecule (£ « 34 meV) and 
the crystal field, if the lattice distortion is applied. 
The matrix of screened Coulomb interactions in the 



-Kg band has been computed in the two steps. First, 
we derive the interaction parameters between atomic p- 
orbitals, using the constraint-LDA method. It enable 
us to obtain the following parameters of intraatomic 
Coulomb interaction u « 11.37 eV, interatomic in- 
tramolecular Coulomb interaction v ~ 2.52 eV, and in- 
traatomic exchange interaction j w 2.30 eV. Then, we 

TABLE I: Transfer integrals and parameters of Heisenberg 
model (in meV) for undistorted bet phase. The nonequiv- 
alent vectors separating two O2 molecules are denoted as 
bi = (0,a,0), b 2 = (f,f,f), b 3 = (a,a,0), and b 4 =(0,0,c). 
Other matrix elements t"j m are obtained using symmetry 
operations and permutability of orbital indices, t}j—t 2 j. 
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take into account the screening of Coulomb interactions 
between p-orbitals in the Tr g band by other bands, which 
are constructed from the same p-orbitals (i.e., the 
band, which is the bonding combination of p x - and p y - 
orbitals, as well as the a* and c* bands, which have a 
strong weight of p z -orbitals). This part is done in the 
random-phase approximation, which starts with interac- 
tion parameters obtained in constraint-LDA. Then, we 
are able to derive the interaction parameters between 
molecular p-orbitals in the ir g band: the intraorbital 
Coulomb interaction U ~ 3.66 eV and the exchange in- 
teraction J w 0.62 eV. The interorbital Coulomb interac- 
tion U' is related with U and J by the identity U'=U—2J, 
These parameters define the whole matrix [/=||£/ Q /3 7 <5|| of 
screened Coulomb interactions in the Tr g band. 



Since any of i™" 



is smaller than U by at least one or- 
der of magnitude, all transfer integrals can be treated as 
a perturbation, starting with isolated molecular orbitals. 
Then, it is convenient to use the hole representation. The 
advantage is that there is only one hole state associated 
with each molecular site. Therefore, if molecules were 
fully isolated, the holes would not interact with each 
other, and the ground state would be a single Slater de- 
terminant. By denoting the hole-orbitals associated with 
i and j as on and otj , and constructing the two- hole deter- 
minant G(ai,aj)=-^={ai(l)aj(2)—aj(l)ai(2)}, one can 
find the energy gain caused by virtual hoppings in each 
bond [13, El: 



J 



T(ati,aj) = — ( G{ai,ctj) 



V j \jM)(jM\P ] 
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G(ati,aj) 
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where Ejm and \jM) stand for eigenvalues and eigen- 
states of excited two- hole configurations at the site j, 
and Vj is a projector operator, which enforces the Pauli 



principle and prevent any hoppings onto otj [Tjj . Eq. @ 
incorporates some intramolecular correlations beyond the 
Hartree-Fock approximation, which are included into the 
multiplct structure of the excited states. 

First, we discuss the properties of the bet phase of KO2 
in an intermediate-temperature region, where the crystal 
distortion is small and fc^T-C^ (note that £/fc,B~ 395 



which is solved self-consistently for each T and the ex- 
ternal field b, acting on f z . The expression for the 

Curie temperature, (kBTc)~ 1 =J2q(Jo~Jq)~ 1 i can be 
derived by assuming &=(f 2 }=0 or from the divergence of 
Xm=4:fi%(T z ) /b. The spectrum of SW excitations itself is 
rather remarkable (Fig. [5]). All excitations are gapped. 
The gap in the Brillouin zone center is proportional to 
£, which is large on the temperature scale. However, 
all low-energy excitations correspond to the zone bound- 
ary, where w q is largely reduced by AFM interactions 



K), so that the splitting of the molecular Tr g levels is J± Thus, the FM order is expected only at T c ~ 66 K, 



largely controlled by the SOI, which provides the natu- 
ral basis for the hole orbitals {cti}. In this case, each 
\oii) can be presented as a linear combination of only 
two orbitals, |p+f) and which are split off by 

the SOI, and where \p±)—^^=(\p x )±i\p y )) are the com- 
plex harmonics. Now, the spin and orbital variables are 
rigidly bound, and the problem can be reformulated in 
terms of the pseudospin states: \t^) = ^(\p+^)±\p-{)), 

\T±)=±{\ P+ ])±i\p-l)), It+HM), and |r-)=b_|), 
which allow us to make a formal mapping of energies ([2|) 
onto the anisotropic Heisenberg model for the bet lattice 
and the pseudospin 1/2: 



H = - 

where f x , 
trices in 
bitals. 
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and f z are the 2x2 Pauli ma- 
basis of |p+T) and \p-i) or- 
The parameters 



'fi- 
the 



calculated as Uh= 
and 



T( T + T -\_q-( T + T +\ = T( T + rT )_q-( T + t +\ 

2jfj— r T{r^ z ,rJ z )— r T{T^ z ,T^ z ), are summarized in Tabled] 
The direct exchange interactions between O2 molecules 
are considerably weaker and can be neglected [f|. The 
interlayer coupling J" operating between neighboring 
molecules separated by b2 stabilizes the FM ground 
state. Hence, the pseudospin moments are parallel to 
the z axis. Other interactions are antiferromagnetic and 
frustrated on the bet lattice. The expectation values 
of the magnetic moments are given by the following 
identities: (r± \L x +2S x \t±)=(t± \L v +2S y \T y t )=0 and 
(rf- \L z +2S z \T z t )—±2: i.e., the finite magnetization is 
allowed only along the z-axis, while it is totally quenched 
in the perpendicular (x and y) directions. The same 
matrix elements define the g-tensor (e.g., in Xm)- 

Then, the magnetic properties can be easily calculated 
using renormalized spin- wave (SW) theory [19j. The idea 
is to extend the conventional formalism to finite tem- 
peratures by replacing in the expression for the SW fre- 
quency the value of local magnetic moment (in our case, 
the pseudospin moment) by its thermal average (f z ): 
w (l =2{f z ){4-J^)+2b 1 where 4'=Ej J ij and J q is the 
Fourier transform of {Jy}. On the other hand, (f z ) is 
expressed through the averaged number of excited spin- 
waves, yielding the equation (f 2 ) cothwk/2fcBT = 1, 



which is largely reduced in comparison with £/fcs (and, 
apparently, can be further reduced by small lattice dis- 
tortions, which are seen in K0 2 HH). The renormalizes 
SW theory also predicts some small deviation from the 
Curie- Weiss law in the paramagnetic region, and the ef- 
fective moment {n e fi), derived from the slope of will 
generally depend on T. Indeed, /i e ff varies from 2 /is at 
large T till 2.17 /j,b near Tc- This may help to explain 
some deviations of fi e g from its nominal value . 
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FIG. 2: The unrenormalized spin-wave dispersion of 
anisotropic Heisenberg model: lj„=2(Jq — Jq)- The inset 
shows the inverse magnetic susceptibility in the paramagnetic 
region, as obtained in the renormalized spin-wave theory. 

In the high-temperature region, fcgT^, the thermal 
fluctuations will eventually destroy the coupling between 
spin and orbital degrees of freedom, which now can be 
treated as independent variables. Then, the spin Hamil- 
tonian can be derived by averaging the pair interactions 
@ over orbital variables (hJ. Practically, for each pro- 
jection of spin, one can consider the whole family of pos- 
sible hole-orbitals, |aj"'')=cos0j|p x t> D+sin^e"^ \p y ], |) 
(where 0<#i<7r and O<0i<27r), and to average (numeri- 
cally) each T(al' l ,a]' 1 ) over all combinations of (9i,4>i) 



and 



'31 



I. The simplest approximation, which is jus- 



tified for large T, is to assume that all combinations are 
equivalent and to neglect the spacial correlations between 
orbitals. Obviously, the spin system will be isotropic, 
and the parameters of Heisenberg model can be found as 
2Jij—'T^—T^, where and are the averaged val- 
ues of T(aJ, cm and T(aJ, aj), respectively. All param- 
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eters are antiferromagnetic (Table HJ, which naturally ex- 
plain the sign of Tew m the high-temperature region Q . 
The absolute value of Tqw can be estimated from the 
behavior of spin susceptibility (xs)i m the framework of 
renormalized SW theory. Note that \s alone should pro- 
vide a good estimate for the total \ m in the perpendicular 
xy-plane, where the orbital magnetization is quenched. 
This yields Tcw~ —100 K, which was derived from the 
linear interpolation of x^ 1 a t large T . 

Thus, the existence of several paramagnetic segments 
with different Tcw's can be related with a gradual in- 
crease of on-site correlations between spin and orbital 
degrees of freedom, driven by the SOI in the bet phase. 
In the high-temperature region, the spin and orbital vari- 
able are totally decoupled, that explains the AFM char- 
acter of intermolecular interactions. On the other hand, 
the strong spin-orbital correlations in the intermediate- 
temperature regime would lead to a ferromagnetism. 

The crystal distortion plays a sizable role only at very 
low T, as it is clearly manifested in the quenching of 
orbital moments and the observed AFM order Q. Al- 
though exact details of the low-temperature structure 
are unknown, a number of experimental data speak in 
favor of uniform rotations of O2 molecules around one of 
tetragonal axes as the main distortion. Certainly, such 
a distortion will lift the degeneracy of molecular p x and 
p y levels, and lead to a FM orbital order, which typi cally 
coexists with AFM correlations between the spins [1 01 ] . 
However, the situation is not so simple. Let us consider 
a distorted structure, which was obtained by the rigid 
clockwise rotation of O2 molecules around the y-axis by 
30° (the experimental value [H, in the fixed tetrag- 
onal frame, and derive parameters of the model Hamil- 
tonian (p]) for the (deformed) Tr g band. The obtained 
"crystal-field splitting" of molecular ir g levels is about 290 
meV, which exceeds £ by nearly one order of magnitude 
and stabilizes the p^-orbital (opposite to the assumption 
made in the phenomenological theory of magnetogyra- 
tion [HI]). Generally, the transfer integrals are specified 
by all three parameters (t 11 , t 12 , t 22 ) 20]. For the neigh- 
boring 2 molecules, they are (-20,0,74), (-25,0,2), 
(37, -29, -67), and (-43, 35, -21) meV, for b x =(0, a, 0), 
bMa,0,0), b 2 =(f,f,§), and !>$=(-§,§,§), respec- 
tively. Thus, the transfer integrals are strongly 
anisotropic, both in and between the tetragonal xy- 
planes. Since the SOI is no longer important, the energies 
^ can be mapped onto rotationally invariant (J^=j]L) 
Heisenberg model. The natural choice of the hole or- 
bitals, |p y t) or \Pyl), is dictated by the crystal-field 
splitting. This yields the following interaction param- 
eters Jij= -0.77, 0, -0.59, 4, -0.09, and -0.01 meV 
between O2 molecules separated by the vectors b 1; b^, 
and b 4 , respectively. They form quasi- two- 



netic order 2l|. Moreover, the AFM interactions in 
each plane are frustrated. Therefore, the Neel temper- 
ature (Tq) is expected to be small. Indeed, from the 
mean-field analysis of the Heisenberg model we conclude 
that the ground state is a (nearly AFM) spin-spiral with 
Qw (0, MI, 2^1). Then, the renormalized SW theory 
yields Tq« 11 K, being in good agreement with the ex- 
periment [8(. The first-order transition to the FM bet 
phase in a magnetic field is also anticipated. 

In summary, KO2 is the new molecular analog of cor- 
related spin-orbital systems. The geometry of molecular 
orbitals adds new functionalities into this classical prob- 
lem. There is also a number of open questions. Par- 
ticularly, any quantitative theory of magnetogyration is 
missing. Why do the O2 molecules tend to rotate at low 
T? Why do their orientation is different in K0 2 and 
Na02? These problems prompt further research. 

I wish to thank P. Mahadevan for drawing my attention 
to K0 2 , as well as M. Boero and Y. Motome for useful 
discussions. This work is partly supported by Grant-in- 
Aids for Scientific Research in Priority Area "Anoma- 
lous Quantum Materials" from the Ministry of Educa- 
tion, Culture, Sports, Science and Technology of Japan. 



J 2- 



ba, b 3 , 



dimensional networks, where interactions in the planes 
perpendicular to [h 1 x b 2 ] are much stronger than those 
between the planes. This suppresses any long-range mag- 
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